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Abstract 

The rapid increase of the e + e~ — > W + W~ cross section in the threshold 
region provides a method for measuring M\y at LEP2. The dependence of the 
theoretical cross section — including the effects of the finite W width, Coulomb 
interactions and initial state radiation — on the collider energy and on the mass 
and width of the W is investigated in order to deduce the optimal collider energy 
which maximizes the sensitivity to M\y. Statistical and systematic uncertainties 
are studied in detail. The background contribution from QCD four-jet produc- 
tion is estimated. It is shown that by running LEP2 at yfs ~ 161 GeV and 
obtaining a luminosity of order 100 pb _1 per experiment, an uncertainty on the 
measured value of M\y of ±70 MeV should be achievable. 



1 Introduction 



A precision measurement of the mass of the W boson is one of the most important 
physics goals of the LEP2 e + e~ collider. The impact of such a measurement on tests 
of the Standard Model and constraints on new physics beyond is well documented, see 
for example Sec. 10 of Ref. [p]] or Sec. 8 of Ref. [@. 
The current 'world average' value is || 

M w = 80.23 ±0.18 GeV , (1) 

obtained by combining the results from the three pp collider experiments: CDF, DO 
and UA2. The precision will improve with increasing luminosity at the Tevatron 
collider. In fact it is likely that the error will decrease to 0(100) MeV in the next few 
years, and may ultimately be as small as 50 MeV 0. Clearly the aim of LEP2 should 
be to at least match this level of precision. 

At least three methods of measuring Mw at LEP2 have been proposed, see for 
example Refs. [0], [| . The most direct method, and at first sight the most promising, is 
to reconstruct Mw directly from its hadronic decay products using the decay channels 

W + W -> qqqq, (2) 
W+W~ -> qqiv . (3) 

If, for example, an integrated luminosity of / C ~ 500 pb -1 can be obtained at a 
collider energy in the range 170—200 GeV, then the statistical uncertainty for this 
method (combining four experiments) is well below 50 MeV H. Notwithstanding the 
difficulties of reconstructing and measuring hadronic jets, the systematic errors appear 
to be under control as well. There is, however, one important caveat. It has recently 
been argued |J that 'colour reconnection' effects — strong interactions between the 
separate hadronic final states of each W in (0) — may distort the mass measurement. 
The studies which yield small systematic errors on Mw from reconstructing the four- 
jet final state assume that the W + W~ system can be treated as a superposition of two 
separate W decays, an approximation which may not be valid at the necessary level 
of precision. 

A second method of measuring Mw is to use the fact that the end-points of the 
lepton spectrum in W + W~ — > qqiv depends quite sensitively on the W mass: 

E.<E t <E+, E ± = ^ (l ± yfl - AM w /s^ . (4) 

The problem with this method is that the statistical error is large for the total lumi- 
nosity envisaged at LEP2: the net branching ratio for the qqiv final state (~ 8/27, 
for I = e, /i) is quite small, and only a fraction of events close to the end-points are 
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sensitive to M w . In Ref. [Q a statistical error on M w of order 300 MeV was obtained 
using this method. 

The third method, the subject of the present work, is to use the rapid increase of 
the W + W~ cross section at yfs ~ 2Mw to measure the mass. Although the rise of 
the Born (on-shell W) cross section is quite sharp, a ~ /3 where (3 2 = 1 — 4M^/s, 
various effects, for example the finite W width, smear out the threshold cross section, 
and so in practice the procedure is to fit a theoretical expression for the cross section, 
with M w a free parameter, to the data. The advantage of this method is that it is a 
direct 'kinematic' measurement of the W mass, in the sense that the dependence of the 
cross section on Mw near threshold comes mainly from the opening up of the phase 
space rather than from the matrix element. There are, however, several important 
questions to address: (i) is it necessary to run at several collider energes to scan the 
threshold region, or is it better to run at a single 'optimal' energy, (ii) what is the 
resulting statistical precision on the mass measurement, (iii) how accurately must the 
cross section be measured, and how precisely does the theoretical cross section have to 
be known, to keep the systematic uncertainties well below the statistical uncertainty? 

In this paper we attempt to answer these questions. We first of all discuss in some 
detail the status of the theoretical calculation of the e + e" — > W + W~ cross section 
in the threshold region. We derive and use an expression which, we believe, contains 
the bulk of the higher-order effects. We then study the variation of this cross section 
with \/s, Mw, TV, ... to estimate the optimal collider energies, and the corresponding 
statistical and systematic uncertainties. 

There have been several other recent studies along broadly similar lines. In the 
theoretical studies of Refs. || 0, the emphasis was on estimating the uncertainty on 
My/ as a function of the collider energy. More recently, ongoing studies within the LEP 
collaborations have also been reported |9], |10j- These address additional issues such 
as the problem of background contributions to the cross section. Where these studies 
overlap with the present work, there is a large measure of agreement. Earlier work 



11, [T2, 13, 14, 15 1 focused mainly on theoretical issues concerning the calculation of 



the W + W~ cross section at LEP2 energies, and also of the impact on the cross section 
of new physics beyond the Standard Model. 

The present work is closest in spirit to Ref. ||, in which many of the concepts to be 
discussed below were first presented. However our study extends the analysis of Ref. 
in several important ways. We include, for example, the leading Coulomb corrections 
which have a significant effect on the threshold cross section, and consequently on the 
extracted value of Mw- We also attempt to quantify the most important background 
cross section, and assess its contribution to the systematic uncertainty on the W mass. 
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2 Calculation of the W + W cross section 

It is not our intention to perform a complete calculation of the theoretical cross sec- 
tion including, for example, all next-to-leading order electroweak corrections or all 
background contributions. In fact the next-to-leading order electroweak corrections 
for four fermion (i.e. off-shell W + W~) production are only partially known at present, 



see for example Ref. [plj for a detailed and up-to-date review. Our aim instead is to 
use a theoretical cross section which includes the most important higher-order effects 
to quantify the dependence on M^. We will then be able to estimate the effect on the 
extracted My/ value of subsequent small variations in the theoretical prediction. 
We begin by writing the cross section, schematically!] as 

CT = 0"ww + Cbkd j 

o"ww = cr^ W (1 + Sew + <5qcd) , (5) 
where the various terms correspond to 

(i) crj vw : the Born contribution from the 3 leading-order diagrams for e + e~ — > 
W + W~ involving t-channel v exchange and s-channel 7 and Z exchange, calcu- 
lated using off-shell W propagators. 

(ii) <5ew : higher-order electroweak radiative corrections, including loop corrections, 
real photon emission, etc. 

(hi) Sqqb'- higher-order QCD corrections to W + W~ final states containing qq pairs. 
Since we are only concerned with the total cross section in this study, the effect 
of these is to generate small corrections to the hadronic branching ratios which 
are entirely straightforward to calculate and take into account. 

(iv) o"bkd: 'background' contributions, for example from non-resonant diagrams (e.g. 
e + e~ — > /i + z/ M iy~) and QCD contributions e + e~ — > qqgg, qqqq to the four-jet 
final state. 

We will begin by studying (i) and (ii) in detail to identify the dominant effects, and 
later estimate the size of the background contributions. Note that all the various 
theoretical components which we need for this study already exist in the literature. 



1 We note that this decomposition of the cross section into 'signal' and 'background' contributions 
is practical rather than theoretically rigorous, since neither contribution is separately gauge invariant 
nor experimentally distinguishable in general. 
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2.1 The W + W off-shell cross section 

The leading-order cross section for off-shell W + W~ production was first presented in 



Ref. |7| : 



where 



a(s) = dsi ds 2 p{s 1 )p(s 2 ) (To(s,s 1 ,s 2 ) , (6) 

JO JO 

C ^ = lliL t (7) 

P{S) ~nM w (s-M^y + s^ w /M^- {i> 

Note that this expression uses a s-dependent W width, 

T w(s) = Tj2- rv , (8) 

where rv = T W (M^), as advocated in Ref. || for example. We will discuss this 
further below. 

The cross section o"o can be written in terms of the v, 7 and Z exchange contribu- 
tions and their interferences: 

cr (s, si, s 2 ) = — — g [a 77 + a Z z + a~,z + «W + a uZ + a vy ] , (9) 

ZoOTTS S\S 2 

where g 4 = e 4 / sin 4 6y/. For completeness, the expressions for the various contributions 
are listed in the Appendix. The stable (on-shell) W + W~ cross section is simply 

a on (s) = a (s,M^,M^) . (10) 

It will be an important factor in the discussion which follows that near threshold the 
cross section is completely dominated by the t-channel neutrino exchange diagram. 
This leads to an S"-wave threshold behaviour at ~ ft, whereas the s-channel vector 
boson exchange diagrams give the characteristic P-wave behaviour a s ~ /3 3 .Q This is 
illustrated in Fig. 1, which shows the on-shell cross section decomposed in terms of 
the six contributions of Eq. ([]). 

2.2 Higher-order electroweak corrections 

The complete set of 0(a) next-to-leading order corrections to W + W~ production 



has been calculated by several groups, for the on-shell case only, see Ref. ]T6| and 



references therein. There has been some progress [18], |19[ with the off-shell (i.e. four 
fermion production) corrections but the calculation is not yet complete. However 
using the on-shell calculations as a guide, it is already possible to predict some of the 



2 We only consider unpolarized e + e collisions in this study. 
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largest effects. For example, it has been shown that close to threshold the dominant 
contribution comes from the Coulomb correction, i.e. the long-range electromagnetic 
interaction between almost stationary heavy particles. Also important is the emission 
of photons collinear with the initial state e ('initial state radiation') which gives rise to 
logarithmic corrections ~ ah\(s/m 2 e ). These leading logarithms can be resummed to 
all orders, and incorporated for example using a 'structure function' formalism. In this 
case, the generalization from on-shell to off-shell W's appears to be straightforward. 
For the Coulomb corrections, however, one has to be much more careful, since in this 
case the inclusion of the finite W decay width has a dramatic effect. Finally, one 
can incorporate certain important higher-order fermion and boson loop corrections 
by a judicious choice of electroweak coupling constant. Each of these effects will be 
discussed in turn below. 

In summary, certain 0(a) corrections are already known to be large, because their 
coefficients involve large factors like log(s/mg), <J M w /T w , mf/M^r, etc. Once these 
are taken into account, one may hope that the remaining corrections are no larger 
than ax constant. When estimating the theoretical systematic uncertainty on the W 
mass in Section |5], we shall therefore assume an overall uncertainty on the cross section 
of ±1% from the as yet uncalculated 0(a) corrections. 



2.2.1 Coulomb corrections 



This section is essentially a summary of the results of Refs. f20|, |21f, where the Coulomb 
corrections for on-shell and off-shell W + W~ production have been discussed in detail, 
and where a complete set of references to earlier studies can be found. 

The result for on-shell W + W~ production is well known — the correction diverges 
as l/v as the relative velocity v of the W bosons approaches zero at threshold. 
Explicitly, to 0(a), 

cr c (s)=a (s)[l + 5 c (s)}, (11) 

where 

Sc(s) = - 7T , (12) 
^0 

with the relative velocity given by 



^0 = 2^1- (13) 
As the threshold is approached from above, 5c — > oo and the corrections must be 



resummed to all orders, giving [22 



1 + - + ... — 2mr/ "° , , . (14) 

vq 1 — exp(— 2an/vo) 
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Note that a ~ v near threshold and so the Coulomb-corrected cross section is for- 
mally non- vanishing when = 2M\y- 

For unstable W + W~ production the finite decay width Y w screens the Coulomb 
singularity p0| , so that very close to threshold the perturbative expansion in a/v^ 



is effectively replaced by an expansion in aJMw/^w plf - The net effect is a cor- 



rection which reaches a maximum of approximately +6% in the threshold region [23 



Although this does not appear to be large, we will see below that it changes the thresh- 
old cross section by an amount equivalent to a shift in Mw of order 100 MeV. The 
result for off-shell W + W~ production is (cf. Eqs. ( |TT| , |T2"D ) 



o (s 



dsi 



ds 2 p(si)p(s 2 ) ooO, si, s 2 )[l + S c (s, si, s 2 )} (15) 



where ||21|| p| 



5c(s, si, s 2 ) 



a 



7r — 2 arctan 



\k\ 2 — p 2 
2p Re(K) 




'2s(si + s 2 ) - (si - s 2 ) 



-M w (E + iT w ) 



E 



4M^ 



4M, 



vv 



(16) 



(17) 



(18) 



In Ref. [21] the above 0(a) result is generalized to all orders. However the con- 
tributions from second order and above change the cross section by 1% in the 



threshold region [24| and will be omitted in what follows. 



2.2.2 Initial state radiation 

Another important class of electroweak radiative corrections comes from the emission 
of photons from the incoming e + and e~ . In particular, the emission of virtual and 
soft real photons with energy E < uj gives rise to doubly logarithmic contributions 
~ a\n(s/m 2 ) \h(s/uj 2 ) at each order in perturbation theory. The infra-red (lncu) loga- 
rithms cancel when hard photon contributions are added, and the remaining collinear 
(hi(s/m 2 )) logarithms can be resummed and incorporated in the cross section using a 

SB a) 



'structure function' [051 (see also Refs 



(Jr. 



r s ds' 

I — F(x, s)\ x= x- s t/ s cr(s') 

*' s min & 



where 26 



F(x, 



tx l - L S + H 



(19) 



(20) 



3 Note that in Ref. ]|^] a slightly different version of Eq. (|l^) is derived, but the numerical difference 
in the threshold region is very small. See Ref. [pTJ for a critical discussion. 
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S -- 
H -- 

with 



3 9 n 2 
1 + - t + 

4 32 12 




1 ) t + 



1 ) In a; 



+ 



1 + 3(1 -x) 



X 



ln(l - x) H 



t 2 + 



t - 2a 

71 



In — - 1 



mi 



(21) 



(22) 



The S term comes from soft and virtual photon emission, while the ri term comes 
from hard collinear emission. Contributions up to 0(t 2 ) are included, as well as a 
finite soft-gluon contribution at 0(a). 

The incorporation of leading logarithmic effects in this way has been widely dis- 
cussed in the literature, see for example Ref. JTBJ and references therein. In fact the 
form given above — in which the cross section is factorized into a structure func- 
tion multiplying the Born cross section — is strictly only valid for processes where 
there is no transfer of electric charge between the initial and final state particles, e.g. 
e + e~ — > Z — > ff. Only the double logarithmic terms are in fact 'universal', i.e. pro- 
cess independent. In general, the subleading 0(a) corrections do not factorize in the 
manner of Eq. (PU]). One possible refinement applicable to e + e~ — > W + W~ involves 
the 'current splitting technique' of Ref. p9|, in which a charge flow is artificially intro- 



duced into the t-channel neutrino exchange in order that a gauge invariant initial-state 
radiation contribution can be defined and factored out. In fact, close to the W + W~ 
threshold the corrections obtained using the current splitting and universal ISR meth- 
ods differ only at the 1% level [p5| . This difference is in any case of the same order as 
the neglected final state radiation contributions. 



2.2.3 Improved Born approximation 

In the Standard Model, three parameters are sufficient to parametrize the electroweak 
interactions, and the conventional choice is {a, G M , Mz} since these are the three which 
are measured most accurately. In this case the value of My/ is a prediction of the model. 
Radiative corrections to the expression for Mw in terms of these parameters introduce 
non-trivial dependences on m t and M^, and so a measurement of M\y provides a 
constraint on these masses. However the choice {a, G M , Mz} does not appear to be 
well suited to W + W~ production. The reason is that a variation of the parameter 
Mw, which appears explicitly in the phase space and in the matrix element, has to be 
accompanied by an adjustment of the charged and neutral weak couplings. Beyond 
leading order this is a complicated procedure. 

It has been argued P0| that a more appropriate choice of parameters for LEP2 is 



the set {MwjGfj,, Mz} (the so-called G^-scheme), since in this case the quantity of 
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prime interest is one of the parameters of the model. Using the tree-level relation 
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e 2 / sin 2 6 W 



AV2G,M 2 



w 



(23) 



we see that the dominant t-channel neutrino exchange amplitude, and hence the cor- 
responding contribution to the cross section, depends only on the parameters My/ and 
G,, 



• ti- 



lt has also been shown 30 



that in the G^-scheme there are no large next-to- 
leading order contributions to the cross section which depend on m f , either quadrat- 
ically or logarithmically. One can go further and choose the couplings which appear 
in the other terms in the Born cross section such that all large corrections at next-to- 
leading order are absorbed, see for example Ref. [31|. However in this study we are 
only interested in the threshold cross section, and so we will simply use combinations 
of e 2 and g 2 defined by Eq. (p3|) for the neutral and charged weak couplings which 
appear in the Born cross section, Eq. (|). 

In summary, the most model-independent approach when defining the parameters 
for computing the e + e~ — > W + W~ cross section appears to be the G^-scheme, in 
which My/ appears explicitly as a parameter of the model. Although this makes a 
non-negligible difference when calculating the Born cross section, compared to using 
a and sin 2 9w to define the weak couplings (see Table below) , a full next-to- leading- 
order calculation will remove much of this scheme dependence ||30||. 



parameter 


value 


M z 


91.1888 


M w 


80.23 


r z 


2.4974 


T w 


2.078 


a' 1 


137.0359895 


G^ 


1.16639 x 10~ 5 GeV" 2 


sin 2 6 W = sin 2 0$T ff 


0.2320 


m e 


5.1099906 x 10~ 4 


(he) 2 


3.8937966 x 10 8 pb GeV 2 



Table 1: Parameter values used in the numerical calculations. Masses and widths are 
given in GeV. 
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2.3 Numerical evaluation of the cross section 

Figure 2 shows the e + e~ — > W + W~ cross section at LEP2 energies. The different 
curves correspond to the sequential inclusion of the different effects discussed in the 
previous section. The parameters used in the calculation are listed in Table 1. Note 
that both the initial state radiation and the finite width smear the sharp threshold 
behaviour at y/s = 2Mw of the on-shell cross section. The different contributions are 
quantified in Table 2, which gives the values of the cross section in different approx- 
imations just above threshold (\fs = 161 GeV) and at the standard LEP2 energy of 
=175 GeV. At threshold we see that the effects of the intitial state radiation and 
the finite width are large and comparable in magnitude. 



&WW 


y/s = 161 GeV 


y/s = 175 GeV 


a (on-shell, a) 


3.813 


15.092 


(To (on-shell, G^) 


4.402 


17.425 


a (off-shell with r w (M w ), G„) 


4.747 


15.873 


(T (off-shell with T w (s), G^) 


4.823 


15.882 


. . . + 0(a) Coulomb 


5.122 


16.311 


... + 0(a) Coulomb + ISR 


3.666 


13.620 



Table 2: Decomposition of the theoretical e + e — > W + W cross section (in picobarns) 
as defined and discussed in the text, at two LEP2 collider energies. 

Our primary aim is to investigate the dependence of the cross section on M w . We 
do this by taking the solid curve in Fig. 2 (G^, finite width, initial state radiation, 
Coulomb correction) as our 'baseline' prediction. Figure 3 shows the cross sections for 
Mw = 80.0 and 80.4 GeV, normalized to that for Mw = 80.2 GeV, as a function of the 
e + e~ collider energy y/s. As expected, the sensitivity to Mw is greatest near threshold. 
Note that for the on-shell cross section, with no initial state radiation, the lower curve 
would decrease to at y/s = 160.8 GeV and the upper curve would increase to oo at 
y/s = 160.4 GeV. Figure 3 also shows the statistical uncertainty on the cross section 
(~ l/y/ N ev ) which would be expected for J £ = 100 pb -1 at each of four separate 
collider energies. The fact that this uncertainty increases as ^/s decreases is simply a 
reflection of the decreasing cross section. However, it is immediately evident that the 
threshold region offers the highest statistical sensitivity to M W - The point of optimal 
sensitivity will be discussed in the following section. 
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3 Sensitivity of the cross section to Mw 



In this section we investigate which LEP2 collider energy offers the maximum precision 
on Myy from a cross-section measurement. There are several factors to take into 
consideration. The first concerns the statistical uncertainty. Far above threshold 
the number of W + W~ events is large but this is outweighed by the very weak M w 
dependence. Far below the (nominal) threshold there are too few events to achieve 
any precision. The second concerns the systematic uncertainties. The important non- 
W + W~ background processes do not have the threshold energy dependence of the 
signal, and therefore constitute an increasingly large fraction of the event sample as 
the threshold energy is approached from above. The optimal collider energy must 
therefore take all these competing factors into account. 



3.1 Statistical uncertainty 

If / C is available for measuring the cross section with an efficiency e at a collision 
energy y/s, then the statistical error on M w can be estimated as 



AM stat = 



da 



dM 



Aa 



da 



dM 



a 



SC. 



= K 



(24) 



where 



K = 



da 



dM 



(25) 



is a measure of the sensitivity to Mw- The dependence of K on y/s is shown in Fig. 5. 
There is evidently a fairly sharp minimum just (~ 500 MeV) above the nominal 
threshold of y/s = 2M W , at which point 



AM 



stat 



90 MeV 



100 pb _ 



(26) 



3.2 Systematic uncertainties 



If the cross section is subject to an overall multiplicative factor C with an uncertainty 
AC (high-order corrections, luminosity, efficiency, ...) then the associated systematic 
error is 



AM sys = 



da 



dM 



a 



AC 



= J 



AC 



where 



J 



da 
dM 



(27) 



(28) 
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The quantity J is shown in Fig. 6 as a function of ^/s. Again we see that the curve 
has a minimum in the threshold region. At the energy (-^/s ~ 161 GeV) at which K 
is minimized we see that J ~ 1.7 GeV, and so 



AM sys ~ 



17MeV 



AC 



C 



x 100% 



(29) 



For an additive uncertainty (for example from the subtraction of a non-VTVT back- 
ground cross section, see below) the contribution to the systematic error is simply 



AM 



sys 



da 



dM 



Aa 



L Aa 



(30) 



The quantity L is shown in Fig. 7 as a function of ^/s in the threshold region. The 
curve has a similar shape to the statistical sensitivity function shown in Fig. 5, although 
the minimum is not quite as sharp. At the energy {^fs ~ 161 GeV) at which K is 
minimized we see that 

" Aa 



AM 



sys 



470 MeV 



1 pb 



(31) 



Finally, there is in principle a systematic uncertainty in M\y arising from an un- 
certainty in the beam energy. It is already apparent from Fig. 4 that a shift in the 
beam energy is equivalent to a shift in M w . Quantitatively, 



AM 



SYS 



da/dE] 



beam 



da/dM 



AE] 



beam 



(32) 



The ratio of derivatives is almost independent of energy in the threshold region and 
has a numerical value very close to 1, i.e. 



AM 



sys 



1.0 AE heam 



(33) 



Current estimates suggest that a Ai? beam of less than 15 MeV can be achieved P2[ , and 
so the systematic error from this source is not expected to be particularly significant. 



3.3 Dependence of the cross section on IV 

In the Standard Model, the decay width of the W is 

T w = [3 + 6{1 + a s (M w )/n + ...}] T(W - 

1+5 su 



G M 3 

[0.1084 ±0.0002] 11 w 



V2 6tt 



(34) 



where a s (My/) = 0.12 ± 0.01 is used in the leptonic branching ratio. The higher order 
electro weak corrections depend weakly on m t and M#, and for canonical values of 
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these 5 SM = —0.0035 PS|, (see also P7|). The main uncertainty therefore comes 
from the value of Mw itself, 



If we use the current world average value (Eq. (p])) for Mw then Tw = 2.078 ± 
0.014 GeV. In contrast, the existing measurements of the W width from the CDF 
collaboration at the Tevatron pp collider are much less precise: 

T w = 2.064 ± 0.061 (stat) ± 0.059(sys) GeV (36) 



from the ratio of W and Z cross sections method [33], and 



T w = 2.11 ± 0.28(stat) ± 0.16(sys) GeV (37) 



from the shape of the lepton p? distribution near the end-point |34j| . Both measure- 
ments are evidently in good agreement with the Standard Model prediction. The 
weighted average of the two measurements is Tw = 2.067 ± 0.082 GeV. 

Concerning the impact of Tw on the determination of Mw, there are two points 
of view that one can adopt. In the context of a Standard Model fit, one can simply 
use the result Eq. (|35|) when studying the variation of the cross section with Mw near 
threshold. In this case, the effect on the cross section as Tw varies with Mw is much 
smaller than the variation with Mw directly, so that for all practical purposes Tw can 
be fixed at its nominal value of 2.078 GeV. 

A less model dependent approach is to study the variation of the cross section 
when Tw is varied in its allowed experimental range, independently of Mw- Figure 8 
shows the cross section near threshold for Ty^ = 2.078 ± 0.1 GeV, as a function of yfs 
for fixed Mw = 80.23 GeV. Note that there is a stable point at a/s ~ 162 GeV where 
the cross section is approximately independent of Tw- From Fig. 8 we can derive?] 



AM 



da/dT 



da/dM 



AT ~ 0.16 Ar , (38) 



where the derivatives are evaluated at a/s =161 GeV. Hence the uncertainty on Mw 
from the current measurement of Tw is ±13 MeV. It is expected, however, that the 
precision of the pp collider Tw determination will improve significantly in the next few 
years 



A final comment concerns the use of the running W width Tw(s) (Eq. (||)) in 
our calculations. It is clear from Table 0, where the cross section is computed for 



4 Note that below a/s ~ 162 GeV, increasing the width has the same effect on the cross section as 
decreasing the mass. 
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both fixed and running widths, that this does make a non-negligible difference in the 
threshold region. However the effect can be largely understood as a shift in the W 
mass, defined as the position of the pole in the propagator: 

(s - M 2 ) 2 + (sY/Mf = (1 + T 2 /M 2 ) (a - M 2 ) 2 + (Mr) 2 (39) 

where M = M(l + T 2 /M 2 )^ 1 / 2 « M - (r 2 /2M) + ...» M - 27 MeV. In fact, this 
is simply a question of the definition of the VF" mass, and is not really an additional 
uncertainty The conventional procedure (used for example in the Z mass measurement 
at LEP) is to adopt the 'running width' definition. 



3.4 Background contributions 

Any additional background contributions to the cross section will reduce the sensitivity 
to Mw There are two main effects: (i) the statistical error on o~ww is increased by 
additional non-Vl^ + W /_ contributions to the event sample, and (ii) an uncertainty in 
the subtracted background cross section becomes an uncertainty on ay/w and hence 
on Mjy. Let us assume that the background cross section CTbkd and its uncertainty 
AcTbkd are small compared to aww- It is straightforward to show that the statistical 
uncertainty on M w is increased to 



AM stat — > AM stat • Jl + , (40) 

V e ' °~ww 

where e is the efficiency for detecting the W + W~ signal. The contribution to the 
systematic error from the uncertainty on the background cross section is given by 
Eq. (|5ID , with Act = A<Jbkd- 



A detailed quantitative assessment of all sources of background events is beyond the 
scope of the present work. It is however possible to give some estimates as to the likely 
size of the most important contributions. We can begin by classifying backgrounds into 
two types: irreducible and reducible. By the former we mean a four-fermion final state 
with the same flavour structure as a W + W~ decay. This is dominated by non-resonant 
electroweak contributions to the 4/ final state. For example, there is a non-resonant 
contribution to e + e~ — ► W + W~ —+ e + vdu from e + e~ — > 7*, Z* — > e + uW~ — > e + vdu 
where the W~ is radiated off an outgoing e~ produced at the 7/Z vertex.0 Such non- 
resonant contributions have been studied in detail in Refs. ]39|] and PD| . The actual 
size of the additional contributions to the WW cross section depends on the particular 
fermion channels, the collider energy and on possible final-state cuts. A semi-analytical 
calculation of the background cross section |39] for a range of 4/ final states (non- 
identical /, no e or u e ) shows that in the threshold region CTbkgd^iO/CT^vK *C 1%, and 

5 This corresponds to a 'reindeer' diagram in the language of Ref. [p9|. 
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therefore the effect on the My/ determination is negligible. It is possible that final 
states such as e — v e W + pose more of a problem, because of the t-channel photon 
exchange contribution. In this case one needs a careful choice of cuts to eliminate the 
potentially large background cross section, see for example Ref. ||0| . 



Reducible backgrounds have a final state which does not correspond to any WW — > 
4/ decay channel but which in practice cannot be distinguished. By far the most im- 
portant of these appears to be QCD four-jet production, e + e~ — > 7*, Z* — > qqgg, qqqq, 
which is a background to e + e~ — > W + W~ — ► qq'qq' '. Since qqgg production domi- 
nates the QCD cross section, there is practically no interference between the QCD 
and W + W~ final states. Although the QCD four-jet background does not appear to 
cause problems in reconstructing My/ from four-jet final states well above threshold 
(see for example Ref. 0), it is potentially more serious for the threshold measurement, 
since the W + W~ cross section is much smaller at threshold and the QCD cross section 
increases with decreasing y/s. 

In order to gauge the size of the QCD background contribution, we have calculated^ 
the cross sections for e + e~ — > qqgg, qqqq in the threshold region. To define a finite-cross 
section we must introduce cuts which keep the massless quarks and gluons energetic 
and non-collinear. We do this with a simple JADE-type algorithm | 42fl , i.e. we require 



that (pi + Pj) 2 > UcutS, where i,j = q,g, and y cnt is a dimensionless parameter. Since 
this is a tree-level calculation at the parton level, it is not straightforward to relate 
it to the measured jet cross section. However this uncertainty can be significantly 
reduced by normalizing the parton-level prediction to the four-jet production rate 
already measured at LEP. Figure 9 shows the fraction of four-jet events measured by 
OPAL |43| at yfs = M z as a function of y cut . The curves are the theoretical predictions 
for two different values of the argument of the strong coupling a s (// = y/s, y/s/ 2) 
chosen to give a reasonable spread of agreement with the data. Figure 10 shows the 
corresponding (absolute) four-jet cross sections at y/s = 161 GeV. In addition to 
the jet separation cut, we have required that in each event there are two jet pairs 
each with an invariant mass within 10 GeV of Mw- The lower dashed lines are the 
QCD predictions corresponding to the 'fits' in Fig. 9. The dotted line corresponds to 
W + W~ — > 4 jets with no cuts, i.e. the total cross section multiplied by a branching 
ratio of 0.46 for the qqqq final state. The solid line is the W + W~ prediction with both 
types of cut included. We see that the QCD curves fall more rapidly with increasing 
y cnt than the W + W~ curve, as expected from the different singularity structure of the 
corresponding matrix elements. 

Taking the average of the dashed curves as our estimate of the QCD background, 
we deduce from Fig. 9 that the square root factor in Eq. fl4*0"|) falls slowly from about 
1.04 at 7/ cu t = 0.01 to about 1.01 at y cnt = 0.07. However this improvement at large 
?/ C ut is more than offset by the overall decrease in the efficiency for detecting the signal, 

6 We use the matrix elements from Ref. [Eli. 
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which is 0.93 at y cnt = 0.01 but only 0.20 at y cnt = 0.07. We conclude that a small 
y cnt value is preferred. If we take y cut = 0.01 as representative, we obtain 

aww(^j, no cuts) = 2.20 pb 

^w(4j) = 1.90 pb (e = 0.87) 

<r QC D(4j) = 0.15 ± 0.03 pb. (41) 

To estimate the overall efficiency we must include also the qq£(= e, \i)v final state. 
Let us assume for simplicity that the efficiency for this channel is 87% also. Then the 
total efficiency can be estimated asQ 

e all = 0.87 x [ BR(W + W~ -> qqqq) + BR(W + W~ -> qqtv) } = 0.65 . (42) 

Note that a 65% efficiency leads, via Eq. (p6|) , to an increase in the statistical error 
by a factor 1.24. If we assume that the background is indeed dominantly due to QCD 
four-jet production, then there is a further increase of 3% from the dilution factor in 
Eq. (g0|). 

4 A strategy for scanning 

It is anticipated that a total luminosity of O(500) pb" 1 will be obtained at LEP2. 
There are strong physics arguments (Higgs and other new particle searches, anomalous 
WWV coupling limits, . . . ) for running at the highest possible energies, yfs ~ 175 — 
200 GeV, rather than at threshold. However, given the importance of a precision Mw 
measurement and remembering that there may be a difficulty in controlling the QCD 
'reconnection effects' for the qqqq final state, it is worth giving serious consideration 
to using some of the available luminosity in the threshold region. Realistically, / C — 
50 — 100 pb -1 could be regarded as reasonable, in which case the uncertainty on Mw 
obtained from the threshold cross section is predominantly statistical. One therefore 
has to think carefully about which collider energy (energies) are optimal. It is clear 
from Figs. 3 and 5 that the only sensible strategy is to run just above the nominal 
threshold. Notice from Fig. 4 that with this amount of luminosity there is no question 
of obtaining any useful information on My/ from the shape of the cross section alone, 
and so running at a single collider energy, or narrow range of energies, is sufficient. In 
doing so we are relying on an absolute theoretical prediction and an absolute cross- 
section measurement. 

Of course the value of the collider energy which corresponds to maximum statistical 
sensitivity depends on Mw itself. However, we already knows the W mass well enough 

7 Note that if electrons and muons from W — > tv decay are also included, this number will increase 
slightly. 
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to be able to specify this energy in advance. The sensitivity parameter K of Fig. 5 has 
a formal minimum at about 500 MeV above the nominal threshold of y/s = 2My/, but 
the variation with y/s is only slight in a region ±0(500) GeV on either side of this. 
Given that when LEP2 starts operation the error on My/ from the pp collider could 
be smaller by a factor of two than its present value of ±180 MeV, one will know the 
optimal collider energy accurately enough. In addition, an initial LEP2 run at a much 
higher energy will allow My/ to be determined even more precisely from reconstructing 
the final state. 

The situation would be very different if one had no a priori knowledge of My?. 
Then it would be necessary to devise a method of locating the collider energy of 
maximum sensitivity without wasting luminosity. This point has been addressed in 
detail in Ref. [f^j. It turns out that a 'data driven scanning strategy' (the BES method) 
can be devised. According to this one makes a sequence of measurements at different, 
carefully specified collider energies, continually updating the My/ determination, until 
one eventually arrives at the optimal energy point just above the nominal threshold. 
The technique has been used successfully to measure m T at lower energies fl4"4]|. 



5 Conclusions 

In this paper we have investigated in detail the measurement of My/ from the W + W~ 
cross section near threshold at LEP2. We have estimated the statistical and systematic 
uncertainties and shown that it will be the former which will dominate the error. It is 
therefore important to choose the collider energy carefully to maximise the sensitivity. 
We have shown that the optimal point is approximately 500 MeV above the nominal 
threshold energy, y/s = 2My/- Based on current My/ measurements, this corresponds 
to y/s ~ 161 GeV. The actual value is not crucial within ±500 MeV or so. Given that 
the existing experimental error on My/ is well within this range, there does not appear 
to be any need to 'scan' the threshold region to locate the most sensitive point. Our 
estimate for the statistical uncertainty on M w at y/s = 161 GeV is given in Eq. (p6|). 
If we assume 

(i) four experiments each obtaining / £ of integrated luminosity, 

(ii) an overall efficiency of 65%, corresponding to detecting the qqqq and qq£(— e, n)v 
final states with an efficiency of 87%, 

then our estimate is 

i 

(43) 

Depending on the size of the background cross section — see below — this number 
may increase by a few per cent. 



AM 
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56 MeV 



100 pb~ 
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sour cp 


AMw 


HO corrections 


Acr/af% N ) x 17 MpV 


luminosity 


A£/C{%) x 17 MpV 


pfficipncy 


Ae/e(%) x 17 MpV 


background 


Act x 470 MpV/pb 


beam pnprgy 


1.0 AEbeam 


dpcay width 


0.16 AT W 



Tablp 3: Summary of systpmatic uncprtaintips on M\y from a mpasurpmpnt of thp 
W + W~ cross spction at y/s ~ 161 GcV. 



Wp havp also invpstigatpd various sourcps of systpmatic uncprtainty, and our rpsults 
arp summarizpd in Tablp |3|. There does not appear to be any single dominant source. 
If wp assump 

(i) by thp timp LEP2 comes into operation thp thporptical cross section will bp 
known to 1% accuracy, which is hopefully conscrvativp, 

(ii) an uncprtainty of 1% on both thp pfficipncy and thp luminosity, 

(hi) an uncertainty on thp background cross section of Acibkd = ±0.03 pb, see 
Eq.flD, 

(iv) a beam energy uncertainty of 15 MeV, which again is probably conservative, 

(iv) either the W width is calculated in the Standard Model or will be measured with 
equivalent precision before LEP2, 

then from Table § our estimate for the overall systematic uncertainty is0 

AM sys = 36 MeV. (44) 

We cannot of course claim to have done a comprehensive study of all of the above 
effects, statistical or systematic, theoretical or experimental. However we have tried 
to use reasonable estimates to set benchmarks against which future, more complete 
studies can be compared. It should be clear from the previous sections where more 
work needs to be concentrated. Our treatment of background contributions has been 
superficial, to say the least. Nevertheless, our overall conclusion is that the idea of 
running LEP2 at the W + W~ threshold in a dedicated attempt to measure Mw from 

8 combining the individual contributions in quadrature 
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the cross section is worthy of serious consideration. With an overall luminosity of 
order 100 pb _1 per experiment, the target error on Mw of ±50 MeV could well be 
achievable. 

A final comment concerns the impact of new physics beyond the Standard Model on 
this method of measuring the W mass. As already stated, we are relying on an absolute 
cross-section calculation, and therefore any additional new-physics contributions will 
distort the M\y measurement. In particular, we have shown that at \fs =161 GeV 
a small shift in the cross section by an amount Acr new is equivalent to a shift in Mjy 
of 0.47 MeV fbT 1 x Aa new . However we have seen (Fig. 1) that at threshold energies 
the t-channel neutrino-exchange diagram is completely dominant. The typical new- 
physics scenarios which have been considered would be expected to affect only the 
s-channel contributions; for example, anomalous r )W + W~ , ZW + W~ couplings or 
new heavy boson production e + e~ — > Z' — > 4/. The corresponding contributions to 
the cross section would therefore be suppressed by at least one additional power of 
(3. In other words, the threshold region is probably safe from most types of new- 
physics contamination, at least at a level which would cause a significant shift in 
Mw- In any case, it is likely that a threshold energy LEP2 run would take place 
after an initial run at higher energy, in which case the presence of new physics would 
almost certainly have already been detected. Nevertheless, a more detailed study of 
new-physics contributions in the threshold region would be very desirable p5| . 
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Appendix 



The following expressions, when substituted into Eqs. (El) and (EI), give the leading- 



order off-shell e + e — > W + W cross section. The results are taken from Ref. 1117 
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8sin 4 6W Gi(s,si,s 2 ) 
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O-vv = G 2 (s,s 1 ,s 2 ) 

duZ = —(V e + 0>e) 7 
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where a e = 1, v e = 1 — 4 sin 2 6^ and 

A(s, Si, s 2 ) = s 2 + s 2 + s 2 , - 2(ssi + SS 2 + SiS 2 ) 
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Figure Captions 



[1] Decomposition of the Born e + e~ — > W + W~ (on-shell) cross section into the 
contributions from t-channel neutrino and s-channel 7*,Z exchange and their 
interferences. 

[2] The cross section for e + e~ — *> W + W~ in various approximations: (i) Born 
(on-shell) cross section, (ii) Born (off-shell) cross section, (iii) with first-order 
Coulomb corrections, and (iv) with initial-state radiation. The parameter values 
are given in the text. 

[3] The cross sections (off-shell, with Coulomb and ISR corrections) for My/ = 80.0 
and 80.4 GeV, normalized to that for My/ = 80.2 GeV, as a function of the 
e + e~ collider energy ^fs. Also shown are 'data points' with statistical errors 
corresponding to / C — 100 pb -1 at each energy. 

[4] The cross sections (off-shell, with Coulomb and ISR corrections) for My/ = 
78.8 + 0.2n GeV (0 < n < 4) as a function of the e + e~ collider energy y/s in the 
threshold region. 

[5] The statistical sensitivity factor K = \fa \da / 'g?M| _1 as a function of the collider 
energy y/s. The arrow indicates the position of the nominal threshold, yfs = 
2M W . 

[6] The multiplicative systematic sensitivity factor J = a \da / dM\~ l as a function 
of the collider energy \fs. The position of the nominal threshold is indicated. 

[7] The additive systematic sensitivity factor L = \da / dM\~ l as a function of the 
collider energy \fs. The position of the nominal threshold is indicated. 

[8] Dependence of the threshold cross section on the width IV- 

[9] The fraction of four-jet events measured by the OPAL collaboration [43| at 



yfs = Mz as a function of y cut . The curves are the theoretical predictions 
using the 0(a 2 s ) qqqq and qqgg matrix elements with two different values of the 
argument /i of the strong coupling a s . 

[10] Four-jet cross sections from W + W~ and QCD four-parton production at y/s = 
161 GeV, as a function of the jet separation parameter y cnt . An additional 
cut requires that in each event there are two jet pairs each with an invariant 
mass within 10 GeV of My. The lower dashed lines are the QCD predictions 
corresponding to the curves in Fig. 9. The dotted line corresponds to W + W~ — > 
4 jets with no cuts, and the solid line has both types of cut included. 
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